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^SJ ' Abstract. Cosmic ray air showers have been known for over 30 years to emit pulsed radio emission in the frequency 

range from a few to a few hundred MHz, an effect that offers great opportunities for the study of extensive air 
showers with upcoming fully digital "software radio telescopes" such as LOFAR and the enhancement of particle 
^ ' detector arrays such as KASCADE Grande or the Pierre Auger Observatory. However, there are still a lot of open 

questions regarding the strength of the emission as well as the underlying emission mechanism. Accompanying 
the development of a LOFAR prototype station dedicated to the observation of radio emission from extensive air 
showers, LOPES, we therefore take a new approach to modeling the emission process, interpreting it as "coherent 
geosynchrotron emission" from electron-positron pairs gyrating in the earth's magnetic field. We develop our 
' model in a step-by-step procedure incorporating increasingly realistic shower geometries in order to disentangle 

, the coherence effects arising from the different scales present in the air shower structure and assess their influence 

on the spectrum and radial dependence of the emitted radiation. We infer that the air shower "pancake" thickness 
directly limits the frequency range of the emitted radiation, while the radial dependence of the emission is mainly 
governed by the intrinsic beaming cone of the synchrotron radiation and the superposition of the emission over 
Q ' the air shower evolution as a whole. Our model succeeds in reproducing the qualitative trends in the emission 

^ I spectrum and radial dependence that were observed in the past, and is consistent with the absolute level of the 

C/3 . emission within the relatively large systematic errors in the experimental data. 
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1. Introduction tinning technical difRculties involved in the radio mea- 
surements and problems with the interpretation of exper- 

In the mid-1960s, Jelley et al. (1965) discovered that ex- imental data. Ground-based particle detectors and later 

tensive air showers (EAS) initiated by high-energy cosmic on fluorescence techniques were so successful that activi- 

rays produce strongly pulsed radio emission at frequen- ties concerning the radio frequency measurements of EAS 

cies around 40 MHz. The discovery triggered intensive virtually ceased. As a consequence, the research on radio 

research and in the following years a number of exper- emission from EAS froze on a rather basic level: while 

iments established the presence of radio emission from the empirical data gathered by the different experiments 

EAS over a frequency-range from a few to a few hundred largely discrepant, the theoretical models mentioned 

MHz. (For an excellent review of the historical develop- above adopt over-simplified geometries, do not incorpo- 

ments and results we refer the interested reader to Allan j-^te relevant shower-characteristics such as realistic parti- 

1971.) Parallel to the experimental work, a number of au- cle distributions or stay on a rather qualitative level that 

thors worked on the theoretical interpretation of the emis- does not allow direct comparison with concrete experi- 

sion processes (Kahn & Lerche 1966; Lerche 1967; Colgate ments. 

1967; Castagnoh et al. 1969 and Fuji & Nishimura 1969). Today, over 30 years after the initial success of radio 

In the early 1970s, however, general interest started frequency measurements of EAS, the field is about to ex- 

to focus on other, at the time more promising methods perience its renaissance. The availability of powerful digi- 

for air shower and cosmic ray research because of con- tal data processing techniques and the advent of digital 
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radio-interferometers such as LOFAR^ (Low Frequency 
ARray) offer the reaUstic perspective to use radio fre- 
quency measurements of EAS as a very powerful and cost- 
effective tool that complements the established techniques 
very well. LOFAR, initially conceived for purely astro- 
nomical purposes, but offering incredible flexibility with 
its ability to form multiple simultaneous beams as made 
possible by its implementation in software, thereby builds 
a bridge between radio astronomy and particle physics. 

Radio measurements of EAS share the main advan- 
tage of optical fluorescence techniques: They allow a 
very direct view into the development of the air shower 
and therefore yield information that profoundly simpli- 
fies the interpretation of data gained by ground based 
particle detectors. At the same time, however, they are 
not hindered by the need for superb observing condi- 
tions (clear, dark, moonless nights far away from any 
light pollution) that limits the duty cycle of optical flu- 
orescence detectors to typically less than 10 %. For a 
purely radio-triggered array with a low number of anten- 
nas, radio detection of EAS should be feasible for ener- 
gies ^ 10^'' eV. With large arrays such as LOFAR or in 
combination with external triggering by particle detector 
arrays such as KASCADE/KASCADE Grande (Antoni 
et al. 2003) or the Pierre Auger Observatory (The Pierre 
Auger Collaboration 1996), the study of EAS ranging from 
^ 10^^ eV up to ultra-high energies would be possible 
(Falcke & Gorham 2003). 

To investigate and develop the potential of LOFAR for 
radio frequency observations of EAS, we currently develop 
the LOfar Prototype Station LOPES (Horneffer et al. 
2002), which is dedicated to the measurement of EAS. 
Obviously, its experimental realisation has to be accom- 
panied by a thorough theoretical analysis of the underly- 
ing emission mechanism, since past theories have not been 
developed to sufficient depth for application to a concrete 
experiment such as LOPES. 

In this work, we take a new approach to the theory of 
radio emission from EAS, namely the interpretation of the 
emission process as coherent synchrotron emission from 
electron-positron pairs deflected in the earth's magnetic 
field (or shorter: "coherent geosynchrotron emission"), as 
proposed by Falcke & Gorham (2003); see also Huege & 
Falcke (2002). Other than Suprun et al. (2003), who re- 
cently simulated geosynchrotron emission from EAS with 
Monte Carlo techniques, we pursue an analytical approach 
to get a better understanding of the effects governing the 
emission. 

We describe the basis of our approach in some detail 
in Sec. 2 and derive some observationally relevant quan- 
tities in Section 3. Sec. 4 summarises the characteristics 
of the air shower development that are needed for a re- 
alistic modeling of the emission process. In Sec. 5-8 we 
develop our model for the radio emission from EAS step 
by step with increasingly realistic geometries, which helps 
in understanding the coherence effects that play a role in 

^ http://www.lofar.org 



shaping the emission spectrum and spatial distribution. 
After a short discussion of the results in Sec. 9 we con- 
clude our work in Section 10. 

2. The geosynchrotron approach 

Two main emission mechanisms have been proposed in 
the past for radio emission from EAS: Cerenkov radia- 
tion from a charge excess moving with a velocity higher 
than the speed of light in the traversed medium (the 
so-called "Askaryan" mechanism motivated by Askaryan 
1962; Askaryan 1965) and acceleration of charged particles 
in the earth's magnetic field. While the former is dominant 
in case of dense media (Buniy & Ralston 2002; Zas et al. 
1992; Alvarez-Muniz et al. 2000), polarisation measure- 
ments in a number of experiments subsequently supported 
the dominance of the geomagnetic emission mechanism for 
radio emission from EAS in air (e.g., Allan et al. 1969). It 
also seems unavoidable in principle for highly relativistic 
charged particles moving in the earth's magnetic field. 

Coherent geosynchrotron emission from highly rela- 
tivistic electron-positron pairs gyrating in the earth's mag- 
netic field represents an equivalent scenario to that of 
the transverse currents of Kahn & Lerche (1966) (and 
other geomagnetic mechanisms) but is particularly ap- 
pealing because it has the advantage of being based on 
well-studied and well-understood synchrotron theory, an 
excellent starting point for the development of our emis- 
sion model. In the case of radio emission from cosmic ray 
air showers, however, coherence effects as well as non- 
periodic trajectories that are usually not considered for 
synchrotron radiation have to be taken into account. 

In order to assess the coherence effects arising from the 
intrinsic air shower structure, we first analyse the emis- 
sion from a specific point during the air shower evolution, 
namely the point of maximum shower development. Only 
in the last step we integrate over the shower evolution as 
a whole, which is effectively "compressed" into the radio 
pulse that the observer receives since the particles have 
velocities ?; « c. 

At this stage, we do not take into account the 
Askaryan-type Cerenkov radiation. In other words, we set 
the refractive index of the atmosphere to unity. 

2.1. Synchrotron-theory: individual particles 

We base our calculations on the formalism developed in 
Jackson (1975). Any acceleration of a charge gives rise to 
electromagnetic radiation. The emission due to accelera- 
tion in the direction of the instantaneous velocity vector is, 
however, insignificant compared to that caused by the per- 
pendicular acceleration (Jackson 1975). As a consequence, 
any arbitrary particle motion, including the helical motion 
of a charged particle in a homogeneous magnetic field, can 
be approximated as an instantaneous circular trajectory 
with adequate curvature radius. 

Retardation effects caused by the finite speed of light 
give rise to strong beaming effects for highly relativistic 
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Fig. 1. Geometry of single particle synchrotron radiation for 
an observer with line-of-sight vector h enclosing a minimum 
angle Q to the instantaneous particle velocity vector v. The 
equivalent curvature radius is given by p, and the emission can 
be conveniently divided into the components and e.\\. The 
particle trajectory lies in the x-y plane. 



particles. For particles with Lorentz factor 7 the original 
dipole emission pattern is beamed into a narrow emis- 
sion cone of order 7"^ semi-opening angle which sweeps 
over the observer in a very short time, leading to strongly 
pulsed emission dominated by frequency components sig- 
nificantly higher than the particle gyration frequency. 

The geometry of the problem corresponds to Fig. 1 if 
one chooses the origin of the coordinate system to lie in the 
point on the particle trajectory where the angle between 
instantaneous particle velocity vector v and line of sight 
vector n reaches its minimum 0. 

Calculation in the frequency domain circumvents prob- 
lems arising from the retardation effects. Jackson defines 
the quantity uj) as a measure of the frequency com- 

ponent Lo of the electric field normalised to unit solid an- 
gle il. In the far-field limit (distance R to the observer 
large compared to the extent of the particle trajectory, i.e. 
use of Fraunhofer-approximation is possible) A(i2, uj) can 
be approximated and conveniently divided into the two 
perpendicular components and e|| defined in Figure 
1. Retaining the phase information, w) can then be 

written as 



8c7r 



ei(-f-f) [-e\\A\\{uo)±e^A^{uj)\ , (1) 



where the plus-sign is to be used for electrons and the 
minus-sign for positrons, e denoting their unit charge. 
Furthermore 
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where lo = 2ttv denotes the angular frequency correspond- 
ing to the observing frequency v, Ka denotes the modified 
Bessel-function of order a, and the curvature radius of the 
instantaneous circular orbit is given by 



P - 



eB sin a 



(5) 



with magnetic field strength B and pitch angle a between 
the particle trajectory and the magnetic field direction. 

Apart from the adopted far-field approximations, the 
derivation of this result incorporates an integration over a 
highly oscillatory integrand only part of which contributes 
significantly. This integration is usually conducted using 
the so-called "method of steepest descents" also known as 
"method of stationary phase" (Watson 1944). Jackson's 
derivation, although somewhat simplified, is correct as 
long as the observing frequency uj is high compared to 
the gyration frequency of the particles in the magnetic 
field. As the latter is around a few kHz and we are only 
interested in observing frequencies > 10 MHz, the Jackson 
result is well suited as the basis for our calculations. It also 
correctly takes into account that the observer sees only one 
flash of radiation from each particle and not the periodic 
repetition that is associated with synchrotron radiation 
in the classical sense, since the mean free path length of 
the particles of ~ 450 m (at a height of 4 km) is very 
small compared with the length of a full gyration cycle of 
~ 20 km. 

The energy spectrum per unit solid angle of a single 
gyrating particle, correspondingly, is given by (Jackson 
1975) 
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Since the energy spectrum is cx |A(il, w)|2 it grows as 
with particle number N if one assumes fully coherent 
emission. Given a specific distance to the observer R the 
frequency component of the i?-field can be calculated as 



E{R,Lu) = 



1/2 



1 
R 



A{R,u;). 



(7) 



For a given (observer-frame) distribution of gyrating par- 
ticles, the corresponding E{R,u}) can then be superposed 
to calculate the total emission. 

2.2. Synchrotron-theory: electron-positron pairs 

In the air shower, electrons and positrons are created in 
pairs. The symmetry arising from the opposite curvature 
of electron and positron trajectories can lead to a sig- 
nificant simplification of the calculation: For an electron- 
positron pair with perfectly symmetric trajectories with 
regard to the observer, the A^\ contributions from the two 
particles add up to 2A|| , whereas the A± contributions 
completely cancel each other. 
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Fig. 2. Misalignment between the electron and the positron 
in an electron-positron pair no longer allows coherent addition 
of the individual emissions. 



This is, however, an overly special case which does 
not adequately represent the problem we are facing. 
Depending on the direction from which the observer sees 
the particle pair, the cancellation of the A± contributions 
as well as the summation of the A|| contributions arc only 
partial. Furthermore, as the pulses emitted by the rela- 
tivistic particles are very short, there is an inherent coher- 
ence length associated to the emissions of the individual 
particles. If there is considerable misalignment between 
the particles, the resulting phase differences destroy the 
coherence as illustrated in Figure 2. Overall, one would 
therefore have to quantify the coherence losses and in- 
complete summation/cancellation arising from the pairing 
through a form factor. 

A more detailed look at the numbers and characteris- 
tics of the particle distributions in the shower, however, re- 
veals that we may indeed assume full summation and can- 
cellation of A II and A± for an "effective" electron-positron 
pair without introducing too large an error. This approxi- 
mation works well if we no longer look at electron-positron 
pairs that actually form together but rather group pairs of 
positrons and electrons together such that their trajecto- 
ries overlap symmetrically as seen by the observer — and 
if we can accomplish this pairing for the vast majority of 
particles. 

For coherent addition of the positron and electron 
emission to be possible, a significant portion of those parts 
of the particle trajectories from which the observer ac- 
tually receives radiation must overlap. (That part has a 
length of ~ 110 m for 7 = 60, given by the length over 
which the instantaneous velocity vector encloses an angle 
< 7 with the observer's line of sight.) In a typical 10^^ eV 
air shower the shower "pancake", even somewhat before 
and after the shower's maximum development, consists of 
~ 10* particles at any time. Even if the particles were dis- 
tributed homogeneously in the shower pancake, this would 
lead to a particle density of ~ 1000 m^"^. For the realistic 
distributions described in section 4, the densities in the 
dominating centre region are a lot higher. This illustrates 
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Fig. 3. Comparison of (solid) and |^x| (short-dashed) for 
u = 100 MHz, 7 = 60 and B = 0.3 G. Absolute scale is arbi- 
trary. 

that each particle (except in the unimportant outskirts of 
the shower pancake) will a priori have a high number of 
particles in its direct vicinity with which it can be paired. 
The probability that there is a significant overlap between 
the paired particles' trajectories is high because the par- 
ticles' mean free path length of ^ 450 m is considerably 
larger than the aforementioned ^ 110 m of the trajectory 
from which the observer receives radiation. 

Whether a consequent pairing with symmetric trajec- 
tories is possible, however, depends critically on the direc- 
tion distribution of the particles' instantaneous velocity 
vectors. Throughout this work we assume a (5-distribution 
of the particle velocity directions at any given point in 
the shower shell, as we choose the initial velocity vectors 
to point radially away from the spherical shower surface. 
In this situation, the pairing of particles with symmetric 
trajectories becomes simple as long as one allows pairing 
between positrons and electrons from generations of par- 
ticles with a certain net offset in generation time. 

In this scenario of high particle density and S- 
distribution of velocity directions, the emission from an 
"effective" electron-positron pair can therefore be approx- 
imated as that from a pair with perfectly symmetric tra- 
jectories: 

£;p(i?,.). ^e'(-f-t)(-e||)^||(.). (8) 

The fact that |y4||| > \A±\, especially for small 6 where 
most of the radiation is emitted (Fig. 3), furthermore 
demonstrates that A± is not dominating the emission, 
anyway, and therefore gives further confidence in the ap- 
proximation. 

Effectively, this result allows us to drop the differen- 
tiation between positrons and electrons and to consider 
only generic "particles" hereafter. The spectrum emitted 
by such an individual particle then corresponds to: 
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Superposition of these spectra for all particles in the 
shower, correctly taking into account the phase differences 
arising from their relative positions, then yields the emis- 
sion from the air shower as a whole. 



3. Observational quantities 

We present a number of relations of the previous results 
to observational quantities: 

3.1. Pulse reconstruction 

The timc-dcpcndcncc of the electromagnetic pulse corre- 
sponding to a given spectrum E{R,lo) can easily be re- 
constructed for a specific receiver bandwidth by an inverse 
Fourier-transform of the remaining spectrum. Hence, if the 
frequency characteristic of the receiver is given by 
the time-dependence of the electric field E{R,t) can be 
calculated as 



E{R,t) = -^[ b{u)E{R, 



uj) e-''^* doj, 



(10) 



where E{R, —u) is given by the complex conjugate of 
E{R,Lo). 

3.2. Conversion of \E{R,lo)\ to 

In the works of the 1960ies and 1970ies, the strength of 
the measured radio emission was usually denoted with a 
quantity in imits of /iV m^'^ MHz"-'^, which was defined 
as the peak electric field strength during the pulse divided 
by the effective bandwidth of the receiver system Au. In 
practice, the total pulse amplitude (in V) at a given ob- 
serving frequency v was derived from the photographed 
oscilloscope traces of the two polarisation directions and 
then converted to an electric field strength (in V/m) tak- 
ing into account the receiver and antenna gain. This field 
strength, representing the projection of the electric field 
vector on the horizontal plane, was then "back-projected" 
to yield the field strength in the plane perpendicular to 
the shower axis and the magnetic field, in which the elec- 
tric field vector lies (see Eq. (55) for = 0). Division of 
the resulting field strength by the effective bandwidth Av 
of the receiver system then yielded e^. 

To compare our theoretical values of \E{R,uj)\ to the 
experimental results for e^, we analytically reconstruct the 
time-dependence of the electromagnetic field pulse E{t) 
for the case of an idealised rectangle filter of bandwidth 
Ai/, over which \E{R,u!)\ is assumed to be constant, and 
which is centred on the observing frequency v. After time- 
averaging over the high-frequency oscillations. e,y then di- 
rectly follows from the peak field amplitude divided by Au 
and is given by 



/ 1 28 

e^ = \ \E{R,oj)\ « 6.4\E{R,w) 

V TT 



(11) 



3.3. LOPES signal-to-noise 

In order to compare our predictions to the abilities of 
LOPES (or any other generic dipole array), we first esti- 
mate the expected signal-to-noise ratio (SNR) for a receiv- 
ing system consisting of an individual inverted V shape 
dipole antenna with gain G = 1.9 and a receiver incorpo- 
rating a filter with bandwidth Ai/ centred on the observing 
frequency u, a square-law detector (i.e., a detector measur- 
ing the received power) and an integrator which averages 
the signal over a time r. 

The noise level of the receiving system, in our case 
dominated by Galactic noise, can be characterised by the 
noise temperature Tgys ~ Tsky(i^)- Comparison with the 
temperature increase AT corresponding to the power of 
the pulse intercepted by the antenna then yields 



SNR = V2A^ 



AT 



■■sys 



(12) 



where the first factor takes into account the increase of 
the SNR due to the number of independent samples accu- 
mulated in case of bandwidth Au and averaging time r as 
determined by the Nyquist theorem. The energy flux of an 
electromagnetic wave propagating through the (vacuum- 
like) atmosphere is given by the Poynting vector, in SI- 
units and omitting the argument R for the fields being 
defined as 



1 



S{t) = E{t) X H{t) = —E{t) X B{t), 



Mo 



(13) 



where /xq = 47r 10 Vs/Am. ks E LB, it follows that 

1 



\S{t)\ = — \E{tf 



1207rO 



\E{i)\\ 



(14) 



For a point-like source, the effective area of a dipole 
antenna is given by (Rohlfs & Wilson 1996) A^a = 
G}? /A-K = G(? /A-Ku"^, so that it receives the power 



\E{t)Y 



(15) 



where the factor 1/2 is introduced for an antenna measur- 
ing only one polarisation direction of unpolarised radia- 
tion. Averaging over the time r then leads to 



<P>r = 



Gc 



Snu'^fioT 
Gc I" 

StTu'^HOT Jul 



(16) 



where uih/i = 2'k{u ± 1/2 Au) and the last step follows 
from Parseval's theorem as long as the bulk of the pulse 
is sampled in the averaging time r. Assuming that the 
spectrum is fiat over the observing bandwidth Av with a 
value |£?(a;')| = \E{2-ku)\ = const, and using Eq. (11), we 
can write 

<P>T « „ ? \E{2Tru)f 2TrAu 



Sttu'^Hqt 

Gc 

Au^Ijlqt 128"'' 



^Au. 



(17) 
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This averaged power is then directly related to AT via the 
Boltzmann-constant fcs by 



AT = 



<P>T 



knAu ' 

so that from Eq. (12) follows 



(18) 



SNR; 



ttGc 



256 V2 u'^iiokB 




(19) 



sys 



Setting T to a sensible value of 2Ai' ^ , wc get 



SNR 



0.5( — 

1.9 



4000 K 



60 MHz 
35 MHz y V 1 AtV m-i MHz"^ 



(20) 



for a typical LOPES antenna with an observing bandwidth 
of 35 MHz centred on the observing frequency 60 MHz 
and an estimate of Tsky(60 MHz) w 4,000 K (Falcke & 
Gorham 2003). 

For a complete LOPES array consisting of A^ant an- 
tennas, the SNR is then increased by an additional factor 
^l/2Ara„t(A^ant " 1) « ^T72iVa„t for large A^ant- 



4. Extensive air shower properties 

Extensive air showers can be initiated by primary parti- 
cles with strongly differing energy and composition and 
at variable inclination angles. Consequently, their proper- 
tics such as the position of their maximum development 
and the longitudinal and lateral distributions of secondary 
particles can vary significantly. 

Additionally, the simulation of air showers consisting 
of > 10* particles with energies in the MeV range created 
in a cascade initiated by primary particles of energies as 
high as 10^*^ eV, is in itself a very difficult process. There 
are elaborate Monte Carlo simulations such as CORSIKA 
(Heck et al. 1998) and AIRES (Sciutto 1999) which them- 
selves incorporate a number of different models for the 
underlying particle interactions. But although these codes 
are very sophisticated, uncertainties remain, especially at 
the very highest energies that arc out of the reach of ac- 
celerator experiments (see, e.g., Knapp et al. 2003). 

At this stage, however, we do not incorporate the re- 
sults of elaborate air shower simulations. We rather re- 
vert to the widely used analytical parametrisations for 
the longitudinal development and lateral particle distribu- 
tions dating back to Greiscn (1956), Kamata & Nishimura 
(1958) and Greisen (1960) which are admittedly crude, 
but as a first step seem adequate to describe the prop- 
erties relevant to our calculations that an "average" air 
shower would have in case of vertical inclination. (For an 
overview see, e.g., Gaisser 1990.) 



4.1. Longitudinal air stiower development 

The longitudinal air shower development can be 
parametrised by the so-called "shower age" s as a function 
of atmospheric depth X: 



s{X) 



3X/Xo 



3X 



{X/Xo) + 2 ln{EjE„it) X + 2X„ 



(21) 



where Xq =- 36.7 g cm~^ denotes the electron "radiation 
length" in air, ii^crit = 86 MeV corresponds to the thresh- 
old energy where ionisation losses equal radiation losses 
for electrons moving in air, and ATm = Aq ln(_Ep/i?crit) 
marks the theoretical value for the depth of the shower 
maximum in this parametrisation. The shower commences 
at s = 0, builds to its maximum development at ,s = 1 
and then declines over the range s = 1-3. Although orig- 
inally developed for purely electromagnetic showers, the 
formula is suitable to qualitatively describe the average 
development of the "clumpier" hadronic showers as well. 
The theoretical X^ value does then, however, not corre- 
spond to the actual position of the shower maximum. For 
purely electromagnetic showers, the development of the 
total number of charged particles (almost purely electrons 
and positrons) can be described by 



N{s) 



0.31 exp [(A/A,|)(1 - fins)] 



(22) 



as a function of shower age. The predicted value for N 
in the shower maximum (s = 1) is very close to the 
Allan (1971) "rule of thumb" A^^ax = Ep/GeY = 10* 
for a lO^'' eV shower. For the position of the shower maxi- 
mum Amax we refer to the measurements and simulations 
presented in Knapp et al. (2003) that suggest a value of 
^max ~ 630 g cm^^ which corresponds to _Ro ~ 4 km for 
a 10^^ eV air shower and to the works of Pryke (2001) as 
well as Abu-Zayyad et al. (2001). 

4.2. Lateral particle distribution 

The lateral particle density can be described with 
the NKG (Nishimura-Kamata-Greisen) parametrisation, 
which without normalisation corresponds to 



1 



r(4.5 - s) 



M 

r 



27rr(s)r(4.5 - 2s) 



s-2 



r 



s-4.5 



(23) 



To avoid the unphysical singularity of the NKG profile at 
the shower centre we cut off the distribution with a con- 
stant value at radii smaller than 0.1 m. The normalisation 
for the integration is chosen correspondingly (see Sec. 6.3). 

The parameters relevant to the NKG distribution, 
shower age s and Moliere radius tm, show a high de- 
gree of degeneracy. The increase in s during the shower 
propagation broadens the lateral distribution, but at the 
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same time the decrease of tm with increasing atmo- 
spheric density tends to narrow it. One can therefore of- 
ten parametrise a given lateral particle distribution with 
a wide range of different values for s and tm , where tm in 
fact need not be close to the theoretical Moliere radius at 
all (Antoni et al. 2001). We here stick to the theoretically 
motivated values of s = 1 for the shower maximum and 
set tm to the Moliere radius at the atmospheric height 
of the maximum derived from the atmospheric density as 
(Dova et al. 2003) 



^Patm(^o) o « S 

rM{ho) — = 9.6 



Patm(/l) 



Patm(/l) 



(24) 



ever, still trace the distribution of e* correctly because the 
number of e* by far exceeds the number of /i^ . An earlier 
analysis of Volcano Ranch data by Linsley (1986) reflects 
the particle distribution at an altitude of 1,800 m, but is 
based on a very low number of events (especially at the 
radial distances up to a few hundred metres relevant to 
our model) and only determines the shower thickness and 
not the functional form of the arrival time distributions. 
Wc therefore base our model on the Agnetta et al. (1997) 
data and use the Linsley (1986) data only for comparison. 

In Agnetta et al. (1997) the measured arrival time 
distribution at a given radial distance is fitted with a T- 
probability distribution function (F-pdf) defined as 



f{t) ^At"" exp(-Ci), 



(25) 



the form of which (cf. Fig. 4) arises from multiple scatter- 
ing events during the shower propagation. While A only 
comprises a normalisation factor, the fit parameters B and 
C of the F-pdf are directly related to the mean arrival time 
<t> and the standard deviation at of the measured ar- 
rival time distributions (Bury 1975), 



B 



<t> 



1 and C 



<t> 



(26) 
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According to the US standard atmosphere of 1977 as im- 
plemented in CORSIKA (Ulrich 1997) the atmospheric 
density at a height of 4 km corresponds to Patm = 
0.82 mg cm^"^, which in turn yields a Moliere radius of 
TM ~ 117 m. 

4.3. Particle arrival time distribution 

Knowledge of the arrival time distributions of particles in 
the air shower is necessary to parametrise the curvature 
and thickness of the shower front as a function of radial 
distance to the core. Unfortunately, the development of 
the particle arrival time distributions during the shower 
evolution is not well established. Agnetta et al. (1997) 
have analysed Haverah Park data of more than 450,000 
air shower events. These lie in the adequate energy range 
of ~ 10""^^ eV, but were measured at an altitude of 220 m 
and cannot differentiate between e* and /i*. They, how- 



20 ^ r ,30 
t [ns\ 

Fig. 4. F-pdf determining the arrival time distribution of parti- 
cles as measured by Agnetta et al. (1997). Solid: in the shower 
centre, short-dashed: 50 m from shower centre, long-dashed: 
100 m from shower centre. 



The dependence of < i > and at on the radial distance to 
the shower core is modeled by a generalised paraboloid of 
the form 



<t>ir),at{r)=F + G (r/ro) 



H 



(27) 



where tq is set to the Moliere radius at ground level of 
79 m. The parameter sets for <t> (r) and at{r) are listed 
as 



Ft 
Gt 
Ht 



Gey 



(8.039 ± 0.068) ns 
(5.508 ± 0.095) ns 
1.710 ±0.059 



(5.386 ± 0.025) ns 
(5.307 ± 0.032) ns 
1.586 ±0.020. 



(28) 



(29) 



Fitting the arrival time distribution with a F-pdf partially 
cuts off the long tail of particles arriving with very high 
delay. Since the radio emission is, however, dominated by 
the bulk of the particles, the effect is negligible for our 
calculations. 

The thickness of the shower "pancake" is directly de- 
termined by at- The "effective curvature" of the shower 
front is governed by two factors. On the one hand, there 
is a delay of the first particles of the F-pdf arriving at 
distance r from the shower core with respect to the first 
particles arriving in the shower centre. This effect is not 
included in the Agnetta et al. (1997) data. Here we as- 
sume that the delay is negligible for the shower distances 
<^ 100 m that we are interested in (the first particles can, 
with good agreement, be assumed to lie on a flat surface). 
On the other hand, the mean particle delay rises as one 
goes to greater distances from the shower core, a fact rep- 
resented by the increase of <t> (r) . The shower curvature 
determined by <t> (r) can be expressed very well with a 
spherical surface of curvature radius K ^ 2, 300 m, as can 
be seen in Figure 5. 
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Fig. 5. Radial dependence of the particle arrival time distri- 
bution. Solid: shower curvature as given by a spherical sur- 
face with K = 2, 300 m, short-dashed: <t>{r) as given by 
Agnetta et al. (1997), long-dashed: crt{r) as given by Agnetta 
et al. (1997), dash-dotted: at,L{r) as given by Linsley (1986). 



For comparison, we also examine the Linsley (1986) 
parametrisation for at, which is defined as 



CTt,L = Gl (1 + r/rL) 



(30) 



where, for a 10^^ eV vertical air shower, we have Gl = 
1.6 ns, TL = 30 m and JJl = 1-68 ± 0.14. As Linsley 
(1986) does not specify the functional form of the arrival 
time distribution, we assume that it also corresponds to a 
F-pdf. However, since <t> and at are not independent in 
this parametrisation, we have to assume a < t > (r) that 
fits the at (r)-distribution given by Linsley. From the fact 
that 



<t>= VTTb 



(31) 



and the relative constancy of B{r) in the Agnetta et al. 
(1997) data, a natural choice for the distribution is given 

by 



<t>A(r-) 

<t>L{r) = ^ at,L{r), 

<^t.A[r) 



(32) 



where the index A refers to the Agnetta parametrisations 
and L refers to the Linsley parametrisations. 

4.4. Particle energy distribution 

The average energy of the electrons and positrons in an air 
shower corresponds to ~ 30 MeV, i.e. 7 ~ 60 (Allan 1971). 
As a very crude approximation, one can therefore adopt 
a mono-energetic configuration of particles with 7 = 60. 
To illustrate the effects induced by a more realistic parti- 
cle energy distribution, we compare this with a (spatially 
uniform) broken power-law distribution of particle ener- 
gies where d7V/d7 rises linearly with 7, peaks at 70 = 60 
and then declines as 7~^: 



-7 



where we set u — 1, w ~ —2 and 71 — 74.2 which cor- 
responds to a peak of the distribution at 70 — 60. One 
can then calculate the emission of an "energy averaged" 
particle through 



E^{R,uj) =po 



p(7) E{R,Lo) d7. 



where the normalisation constant po is 
1 

Po 



C'^.'^^Pil) d7' 



(34) 



(35) 



This energy integration leaves the number of particles 
unchanged. Note, however, that the total amount of en- 
ergy in the particles varies with changing 7niin and 7niax- 
Additionally, 7min must not be chosen too small as our 
derivations include approximations that are only valid in 
the ultra-relativistic case. In general, the presence of high- 
energy particles amplifies the emission near the shower 
centre, whereas low-energy particles enhance the radiation 
at high distances due to their wider beaming cone. 

We will continue to compare results with both mono- 
energetic and broken power-law particle distributions and 
differentiate the two cases through the absence or presence 
of an additional index 7. Any result calculated for a broken 
power-law distribution, indicated through an index 7, also 
applies to the mono-energetic case if one substitutes the 
energy-averaged Ej{R,u)) by the original E{R,uj). 

5. Coherence: longitudinal effects 

Having established the emission from individual particles 
and the spatial distribution of particles in the air shower, 
we can now calculate the emission from the shower max- 
imum, taking into account the inherent (observer-frame) 
shower structure. The phase differences between the pulses 
from the individual particles lead to strong coherence ef- 
fects that significantly change the spectrum of the received 
emission from that of a fully coherent synchrotron pulse. 
(For a general discussion of coherence effects regarding 
synchrotron radiation see also Aloisio & Blasi 2002.) 

To get a better understanding of the spectral features, 
we start with a strongly simplified configuration that ne- 
glects any lateral structure: We reduce the shower to its 
core. In this approximation, the air shower "pancake" of 
thickness d is collapsed to a one-dimensional line of length 
d. The charges distributed along the line are adopted to 
move synchronously, i.e. the momentum distribution of 
the particles at a given time corresponds to a (5-function. 
The emission from a particle at distance R = fiR from the 
observer is given by E{R,uj) as defined in Eq. (9). For a 
particle offset by a distance x from the line-centre (located 
at Rq ) along the shower core direction I , we therefore have 



E{Ro + xl,ijj) (X 



1 



1 



■xl\ 



exp 



lUj\Rn 



ILO I Rq + xl I 

xh 
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= — exp 
no 



ia;(i?o + x) 



An{uj). 



(36) 



Here we keep the distance of the particle at the constant 
value Rq for the first factor, which only introduces negli- 
gible errors since d <C -Ro- The approximation of Z ~ n in 
the second factor is justified since ^|| only gives significant 
contributions if the directions of I and n enclose angles of 
order 7"^ or smaller. In other words, projection effects 
do not play a significant role because only in the regime 
where they arc very small, we have significant emission 
from the particles. 

Defining the particle distribution function f{x) such 
that f{x) dx = 1, and taking into account the par- 
ticle energy distribution, the integrated emission from N 
particles is then given by 



/+(X 
Nf{x) E^{Ro+xi,uj) 
-00 



dx 



NE^{Ro,iu) / f{x) e" 
NE^{Ro,u}) S{u). 



+ 00 



dx 



(37) 



Note that this basically corresponds to a Fourier trans- 
formation, i.e. the function S{lu) modulating the field 
strength spectrum is given by the Fourier transform of 
the particle distribution function, as in standard diS'rac- 
tion theory. (The energy spectrum is then modulated by 
\Siu;)\\)^ 

We will now compare a number of different distribu- 
tions of particles along the line to better understand the 
coherence effects that arise from longitudinal distributions 
of particles. 



The coherence function then equals 



S{uj) = exp I — 



2 c2 



-U} 



(41) 



i.e. a Gaussian as well, which is clear from the fact that 
the Fourier transform of a Gaussian is again a Gaussian. 

5.3. Asymmetrical T -distribution 

A realistic longitudinal particle distribution is given by an 
arrival-time distribution as specified by Eq. (25) with the 
substitution x = ct, 







exp I 



-C ? 



for X > 
for a; < 



m = 

where from the normalisation of f{x) follows 



c-(i+B)r(i + 5) 



(42) 



(43) 



The corresponding coherence function S{w) is then given 

by 



= (1 + ^ 



4(1+5) 



exp 



i(l + B)arctan(^) 



X exp (— iw <t>A) , 



(44) 



where the last phase factor is needed to "centre" the asym- 
metrical distribution on the curved shower surface to make 
it comparable to the symmetrical uniform and Gaussian 
distribution for the later calculations taking into account 
lateral structure. (The origins of the F-pdfs then again lie 
on a flat surface as discussed in Section 4.3.) 



5.1. Uniform line charge 



5.4. Model calculations 



The easiest case of a line charge is a uniform distribution 
of particles along a line of length d, 



fix) 



1/d 




for 
for 



\x\ < d/2 
\x\ > d/2 



(38) 



Integration over x then leads to the well-known (sinz/z)^ 
modulation of the energy spectrum that corresponds to 
the dififraction pattern of a rectangular opening. 



5 1 



dx = 



sin {d(jj 1 2c) 
rfa;/2c 



(39) 



5.2. Gaussian line charge 

A more realistic case is that of a Gaussian distribution 
of particles along the line. The width of the distribution 
is set by the standard deviation of the Gaussian a (the 
FWHM then corresponds to \/4 ln4 a w 2.35 a), with 
the distribution being defined as 



f{x) = 



1 



fTV27r 



2^ 



(40) 



The results derived so far allow us to perform a number 
of model calculations in order to analyse the effects of 
longitudinal particle distributions on the observed spectra 
as well as the dependence of the emission on the observer's 
radial distance from the shower core. Where no analytical 
result was presented, integrations and other calculations 
are done numerically. We model the maximum of a vertical 
air shower with primary particle energy Ep = 10^^ cV and 
therefore N = 10* particles at a height of i?o =4 km. This 
is a realistic value as outlined in Section 4.1. The earth's 
magnetic field is adopted with a strength of B = 0.3 G 
and, for simplicity and symmetry reasons, assumed to be 
perpendicular to the shower core and thus parallel to the 
earth's surface (a realistic value for Central Europe would 
he B = 0.5 G with declinations around 70°). 

The thickness of the air shower "pancake" is deter- 
mined by the standard deviation at as parametrised in 
Sect. 4.3. To ensure an equivalent width of the F-pdf, the 
uniform line charge and the Gaussian line charge, we nor- 
malise the latter distributions such that they have a stan- 
dard deviation of cat- For the Gaussian distribution, this 
corresponds to a = cat. The uniform line charge must 
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Fig. 6. \E{R, 2-Kv) |-spectrum in the centre of the area illmni- 
nated by the maximum ofalO^^eV air shower with _Ro = 4 km 
and 7 = 60. Solid: full coherence, short-dashed: uniform 5.6 m 
line charge, long-dashed: Gaussian line charge with g — 1.61 m, 
dash-dotted: asymmetrical F-distribution with cat = 1.61 m 



Fig. 7. Radial dependence of \E{R, 2tt 50 MHz)| for the max- 
imum of a lO^'^ eV point source shower with Ro = 4 km for 
the 7 = 60 case (solid) and for broken power-law distributions 
from 7 = 5-120 (short-dashed), 7 = 5-1000 (long-dashed) and 
7 = 5-10000 (dash-dotted). 



be set to a total length of d = 2\/3cat- Evaluated in the 
shower core, ctTt(O) — 1.61 m, which results in c? = 5.6 m. 

Fig. 6 compares the spectral modulations arising from 
the different longitudinal particle distributions. If the par- 
ticles radiated fully coherently — i.e. moved "as one par- 
ticle" on the exact same trajectory — the field strength 
spectrum of the emission would simply be a synchrotron 
spectrum enhanced by a factor TV. The coherence effects 
modulate this spectrum by the coherence function S{lu). 
In the case of the uniform line charge, we see the first co- 
herence minimum at w 54 MHz, which corresponds to c/d. 
The Gaussian line charge spectrum does not exhibit such 
a sharp minimum, but is strongly attenuated at higher fre- 
quencies. The asymmetrical F-pdf lies between these two 
simplified models. 

Obviously the longitudinal effects very strongly modu- 
late the emitted spectrum at high frequencies 50 MHz) 
and therefore are an important limiting factor for the 
choice of a suitable observing frequency. The thickness 
of the air shower "pancake" has a very profound and di- 
rect influence on the emitted radiation and could therefore 
be probed directly through observations of radio emission 
from EAS at frequencies > 50 MHz. 

Another interesting characteristic of the radiation is 
its radial dependence at a given frequency as illustrated 
by Fig. 7 for the case of full coherence without any lon- 
gitudinal distribution, i.e. for particles concentrated in a 
point source. In this case, the associated emission pattern 
is purely governed by the inherent emission pattern of the 
synchrotron pulses. The extent of the illuminated area on 
the ground is an important characteristic that ultimately 
limits the probability to detect very scarce ultra-high en- 
ergy cosmic ray air showers with a given collecting area. 

As expected, adoption of the broken power-law distri- 
bution of particle energies influences the radial emission 
pattern. Going from the mono-energetic case to the broken 
power-law energy distribution mainly amplifies the emis- 
sion in the centre region due to the presence of high energy 



particles that radiate more strongly, but into a smaller 
beaming cone. At the same time, the low-energy parti- 
cles amplify the emission at very high distances due to 
their wider emission pattern. The drop in the number of 
medium-energy particles is correspondingly reflected in a 
drop of the emission on medium scales. Obviously, there is 
only negligible difference when increasing the upper limit 
7max fi'om a value of 1000 to higher values such as 10000. 
For the remaining calculations, we therefore adopt a dis- 
tribution with 7 in the range 5-1000 to minimise compu- 
tation time. 

6. Coherence: lateral effects 

After having analysed coherence effects arising from lon- 
gitudinal distributions of particles, we now take a look at 
the influence of the lateral structure of the air shower on 
the radio emission. This we accomplish by "smearing out" 
the line charge considered so far over a segment of a spher- 
ical surface with (for the moment) constant thickness d. 
Inside this "shell" we continue to consider the types of 
longitudinal particle distributions introduced in Section 
5.4. 

6.1. Geometry 

The geometry of the air shower maximum is deflned as 
in Fig. 8 and characterised by the curvature radius of the 
shower surface K, the shower shell thickness d and the 
shower inclination angle rj. The observer is positioned on 
the X-axis at a minimum distance Rq from the shower sur- 
face, with an inclination angle 7?o to the shower core. The 
magnetic fleld strength B, inclination rjB and azimuthal 
direction ips determine the conflguration of the earth's 
magnetic fleld. 

To derive the total emission from the air shower max- 
imum, we now have to integrate over the shell and hence 
must relate E{R, u) and consequently the quantities go- 
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Fig. 8. The geometry for the air shower maximum. 

ing into E{R,lu) to the position on the surface as given 
by I? and ip. We refer the reader to the appendix for the 
details of these calculations. 



6.2. Approximations 

In order to facilitate the integration over the shower shell 
we apply a number of approximations. First, we sum the 
contributions from the different regions of the air shower 
maximum in a scalar way, i.e. we do not evaluate Eq. (55) 
but set e||(i?, (p) = e||(0, 0), which is justified due to the 
minute changes in the direction of e|| (i9, ip) over the shower 
surface. The general polarisation direction of the radia- 
tion then is perpendicular to both the shower axis and 
the magnetic field direction. Second, as pointed out in the 
appendix, we assume that the instantaneous velocity vec- 
tors of the particles in the shower shell are perpendicular 
to the shower surface at the moment corresponding to the 
origin of Figure 1 . This again corresponds to a (5- function 
for the distribution of the particle momenta, and in this 
strict sense, the minimum angle to the line-of-sight 9{'d, ip) 
is given by Eq. (54). 

Adoption of this 9, however, yields a very conservative 
estimate for the emission. 6{-d, ip) as calculated by Eq. (54) 
overestimates the minimum angle to the line-of-sight as 
defined in Fig. 1 in case of a particle trajectory bending 
towards or away from the observer, where a significantly 
reduced 9 is reached during a later or earlier position on 
the particle trajectory. The amount of "compensation" in 
9 attainable by this effect is considerable since the ratio of 
mean free path length to curvature radius of 7 = 60 elec- 
trons is ^ 450m/3400m« 87"^. 9 is therefore effectively 
reduced to its irreducible component given by 



sin9{-d,ip) = B ■ h{d,ip). 



(45) 



Adoption of this value for 9{{}, ip) yields a more realistic 
estimate of the emission from the air shower shell, and at 
the same time takes into account the asymmetry of the 



emission pattern in ip that arises from the magnetic field 
configuration. Without a more precise criterion for the 
maximum compensation achievable, however, the radial 
dependence of the emission pattern at very high distances 
is obviously no longer valid. We therefore continue to work 
with both the conservative approach using Eq. (54) and 
the "reduced 0" definition in order to compare the two 
cases. 

The change in R associated to the adoption of the "op- 
timum position" on the particle trajectory is negligible 
because of the following reasons: 

— the compensated angles are small, therefore the 
changes in R are small 

— additional attenuation/amplification through the 1/R- 
term is thus negligible 

— there is no significant change of phase since the particle 
velocity v ~ c and the trajectory is only mildly curved. 

6.3. Integration 

Using these geometrical relations, the integrated spectrum 
of the emission from the air shower maximum with N 
particles can be calculated as 



= po dip di) K^aini} 

Jo Jo 

X PNKG(r(i9,¥')) E^{R{i^,ip),uj) 



with the normalisation factor 



po^N 







2ir pr-^i/K 

dip I d-d K'^ siiid p^Y.G{T{d,ip)) 
'0 



(46) 



(47) 



Cutting off the integration dX d = ryijK significantly re- 
duces computation time while giving acceptable accuracy 
as ^ 80 % of the particles are included in this region. The 
remaining particles are redistributed over the integration 
region by the normalisation according to the NKG-profile, 
which might lead to a slight overestimation of the emission 
strength near the shower centre. 

6.4. Model Calculations 

We again examine the frequency and radial dependence of 
the emission to study the effects introduced through the 
lateral particle distribution. The basic parameters adopted 
are the same as in Sec. 5.4, but we assume a broken power- 
law particle energy distribution from 7 — 5-1000 for all 
calculations. The curvature radius of the shell is adopted 
a,s K = 2, 300 m and the Moliere radius set to tm ~ 117 m 
as discussed in Sects. 4.3 and 4.2, correspondingly. 

In Fig. 9 we plot the spectrum received by an observer 
in the centre of the area illuminated by the air shower 
maximum, considering the same set of longitudinal parti- 
cle distributions as before. The spectra look very similar 
to those of a line charge, but are attenuated additionally 
at high frequencies. 
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Fig. 9. |i/(J?, 27r!/) [-spectrum at the centre of the area illu- 
minated by the maximum of a 10^^ eV air shower with re- 
alistic lateral distribution, 7?o = 4 km and a broken power- 
law energy distribution from 7 = 5-1000. Solid: full longi- 
tudinal coherence, short-dashed: uniform 5.6 m longitudinal 
distribution, long-dashed: Gaussian longitudinal distribution 
with G = 1.61 m, dash-dotted: longitudinal F-distribution with 
C(Tt = 1.61 m. For comparison: fully coherent case without lat- 
eral distribution (dotted) 
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Fig. 11. Same as Fig. 10 but with "reduced and magnetic 
field parallel to the direction of the observer. 
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Fig. 10. Radial dependence of \E{It, 2nv) \ for the maximum of 
a 10'^ eV air shower with full longitudinal coherence, realistic 
lateral structure, "conservative 9" approach, Rq = A km and a 
broken power-law energy distribution from 7 = 5-1000. Solid: 
= 50 Mhz, short-dashed: = 75 Mhz, long-dashed: v — 
100 Mhz 



A more interesting result is illustrated by Fig. 10 which 
demonstrates the effect of a purely lateral particle distri- 
bution on the radial dependence of the emission in the 
"conservative 0" scenario, completely ignoring any longi- 
tudinal effects. The lateral structure introduces a mod- 
ulation of the radial dependence, caused by interference 
of emission from opposite ends of the shower "disk" . 
For higher frequencies and, correspondingly, shorter wave- 
lengths, the interference minima move to smaller radial 
distances. 

In comparison. Figs. 11 and 12 show the radial depen- 
dence in case of the "reduced 9" calculations. The inter- 
ference effects are somewhat washed out and the overall 
emission level is higher. As expected, there is a drastic 
asymmetry between the emission pattern in the directions 
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Fig. 12. Same as Fig. 10 but with "reduced 9" and magnetic 
field perpendicular to the direction of the observer. 



parallel and perpendicular to the magnetic field. In case 
of Fig. 12, where the observer is positioned in a direc- 
tion perpendicular to the magnetic field, 9 is basically re- 
duced to zero even for distances of a few hundred metres 
(9 <; 87^^ as explained in Section 6.2). Correspondingly, 
the emission pattern is only very slightly attenuated even 
at high distances. 

7. Flaring disk 

We now combine the results derived so far to obtain a more 
realistic model of the emission from the maximum of an 
extensive air shower: a "flaring" disk. In other words, we 
adopt the same geometry as specified in Sec. 6, but now 
vary the thickness of the disk as a function of position 
(??, if) on the shower surface in the form of the varying 
asymmetric F-pdfs parametrised as in Section 4.3. This 
geometry therefore correctly takes into account the curva- 
ture, the lateral and the longitudinal structure of the air 
shower maximum. 

Fig. 13 again shows the spectrum emitted by the air 
shower maximum as a realistically flaring disk according 
to the Agnetta et al. (1997) and Linsley (1986) parametri- 
sations. As expected, the spectrum emitted by the Linsley 
flaring disk extends to higher frequencies than the one gen- 
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Fig. 13. \E{R, 27ri/)|-spectrum at the centre of the area illu- 
minated by the maximum of a 10^'' eV air shower with flaring 
F-pdf, -Ro = 4 km and a broken power-law energy distribution 
from 7 = 5-1000. Solid: flaring Agnetta et al. (1997) lateral 
distribution, short-dashed: flaring Linsley (1986) lateral distri- 
bution 
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Fig. 14. Radial dependence of \E(R,2-kv)\ for the maximum 
of a 10" eV air shower with flaring Agnetta et al. (1997) F-pdf, 
7?o = 4 km and a broken power-law energy distribution from 
7 = 5-1000. Solid: = 50 Mhz, short-dashed: v = Mhz, 
long-dashed: v = 100 Mhz, upper/lower curves for "reduced 
0" perpendicular/parallel to magnetic field direction 



erated by the Agnetta flaring disk because of the lower 
thickness in the shower centre where most of the particles 
reside (cf. Fig. 5). 

The radial dependence at different frequencies is once 
again shown in Figure 14. Comparison with the cor- 
responding diagrams for the purely lateral distribution 
shown in Figs. 10-12 shows that the overall emission level 
drops as the observing frequency is increased due to the 
dampening of higher frequencies by the longitudinal par- 
ticle distribution. Additionally, one can again observe a 
"smearing out" of the interference minima. As a conse- 
quence, the "conservative 6*" and the "reduced 6*" with 
observer parallel to the magnetic field calculations yield 
almost identical results. 

In Fig. 15 we have reconstructed the pulse generated 
by the flaring Agnetta disk as it would be measured by 
a receiver with a given bandwidth using Eq. (10). The 
pulse amplitude drops noticeably when the observer moves 




Fig. 15. Reconstructed pulses emitted by the maximum of 
a 10" eV shower with flaring Agnetta et al. (1997) F-pdf, 
broken power-law energy distribution from 7 = 5-1000 and 
-Ro = 4 km, using an idealised rectangle filter spanning 40- 
160 MHz and "conservative scenario. Solid: centre of illu- 
minated area, short-dashed: 100 m from centre, dash-dotted: 
250 m from centre 



from the centre of the illuminated area on the ground to 
a distance of 100 m, and is already quite diminished at 
a distance of 250 m, as expected for the "conservative ^" 
approach. The pulse length of ~ 8 ns is a result of the 
filter bandwidth of 120 MHz, i.e. the pulse is bandwidth- 
limited. 



8. Integration over shower evolution 

The last step in modeling the total air shower emission 
is to integrate over the shower evolution as a whole. This 
can be done in a very simplifled fashion by approximat- 
ing the shower evolution with a number of discrete steps. 
The characteristic scale for these steps is given by the 
"radiation length" of the electromagnetic cascades in air, 
Xq = 36.7 g cm~^, corresponding to w 450 m at a height 
of 4 km. One can therefore discretise the shower evolution 
into "slices" of thickness Xq, assuming these contain in- 
dependent generations of particles and therefore radiate 
independently. 

The emission from each of these slices is calculated as 
that from a flaring disk, taking into account changes of s, 
i?o, "i^Oi "^M and N correctly through the relations given 
in Section 4 and reverting to the "conservative 0" defini- 
tion to be able to correctly calculate the emission at great 
angles. Superposition of the individual slice emissions, cor- 
rectly taking into account the phases arising from arrival 
time differences, then leads to the total emission of the 
shower. 

Slices far away from the observer are attenuated both 
due to the high distance and the decreasing number of 
particles N . The concrete number of far-away slices taken 
into account is therefore uncritical. The situation is differ- 
ent for the slices close to the observer. In their case, the 
attenuation through the decreasing number of particles 
ISl is more than compensated by the decreasing distance 
to the observer. In fact, the slices closest to the observer 



14 



T. Huege and H. Falcke: Radio Emission from Cosmic Ray Air Showers: 



> 

3, 



3 



10 



. 1 



.01 









■ 














\ ^ 






\ 






\ 






\ 






\ 






\ 






\ 






\ 






. \ 





10 



100 
u [MHz] 



100 



10 w 



> 

3, 



lOr 



. 

1000 



> 

3, 



0.1 



" 3 0.01 

til 




100 200 300 400 500 
distance from shower centre [m] 



10 



. 1 



600 



> 
3. 



0.01 



Fig. 16. |i?(J2, 27r!/)|-spectrum of a full lO^'^ eV air shower 
with flaring Agnetta et al. (1997) F-pdf, "conservative 6" ap- 
proach, i?o = 4 km and a broken power-law energy distri- 
bution from 7 = 5-1000. Solid: centre of illuminated area, 
short-dashed: 100 m from centre, long-dashed: 250 m from cen- 
tre, black points: rescaled Spencer (1969) data as presented by 
Allan (1971), grey points: rescaled Prah (1971) data 



Fig. 17. Radial dependence of |£:(i?,27r 55 MHz)| for a full 
10'^ eV air shower with flaring Agnetta et al. (1997) F-pdf, 
"conservative 9" approach, Ro = 4 km and a broken power- 
law energy distribution from 7 — 5-1000, data from Allan et al. 
(1970), horizontal lines from top to bottom: emission strength 
needed for a Sa-detection with an individual LOPES antenna 
or an array of 10 or 100 LOPES antennas 



yield the highest contributions of radiation, and the to- 
tal result depends considerably on the number of nearby 
slices taken into account. However, at the same time, the 
illuminated area on the ground, governed by the intrinsic 
beaming cone, becomes very small for the slices very close 
to the observer, especially for the high frequencies where 
the radiation mainly originates from high-energy particles 
with even smaller beaming cones. Except for low frequency 
emission in the centre region of the illuminated area, the 
result for the total emission can therefore be considered 
robust. 

For our vertical 10^^ eV air shower at a height of Rq = 
4 km we add the emission from eight slices above and eight 
slices below the shower maximum to the emission from the 
maximum itself. The closest slice then lies at Rq = 950 m 
from the observer, a distance we do not want to fall below 
because of approximations contained in our calculations 
that are only valid in the far-field. 

The main effect of the integration over the shower evo- 
lution is a boosting of the total emission because of the 
increased total number of particles taken into account, as 
can be seen in the spectra shown in Fig. 16. For frequencies 
of ~ 40 MHz and radial distances of ^ 100 m, the amplifi- 
cation factor corresponds to ^ 10. Apart from the overall 
amplification, the radial dependence is significantly steep- 
ened because the important nearby slices only contribute 
at low radial distances as discussed earlier. This can be 
seen when comparing Fig. 17 with the earlier results for 
the "conservative 6'" case. For illustration purposes, we 
also present a reconstructed pulse from the shower as a 
whole as it would be measured by an observer in the cen- 
tre of the ilhmrinated area using the LOPES bandwidth 
of 35 MHz in Figure 18. 




-40 -20 20 40 60 80 
t [ns] 



Fig. 18. Reconstructed pulse in the centre of the area illumi- 
nated by a full 10^'' eV shower with flaring Agnetta et al. (1997) 
F-pdf, broken power-law energy distribution from 7 — 5-1000 
and _Ro = 4 km, using an idealised rectangle filter spanning 
42.5-77.5 MHz 

9. Discussion 

The calculations presented here represent only a few il- 
lustrative examples of possible configurations of EAS that 
could be calculated with our model. These examples, how- 
ever, already demonstrate the most important dependen- 
cies between shower structure and emission spectrum as 
well as radial emission pattern. 

9.1. Theoretical results 

As expected, the thickness of the air shower pancake, and 
correspondingly in our model the width of the longitudinal 
particle arrival time distributions, is the main factor de- 
termining the position of the high-frequency cut-off in the 
spectrum. Typical longitudinal scales of a few metres lead 
to frequency cut-offs in the 100 MHz regime, which sup- 
ports a choice of observing frequency well below 100 MHz. 
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Due to the strong dependence of the spectral cut-off on the 
shower thickness, radio emission from EAS could be used 
very effectively to probe the longitudinal structure of air 
showers during their evolution, a quantity that is not well 
known at the moment. 

The radial omission pattern is mainly governed by the 
inherent emission pattern of the synchrotron pulses and 
the superposition of the beamed emission from different 
parts of the air shower evolution as a whole. Additionally, 
the lateral extent of the air shower slightly influences the 
size of the illuminated area on the ground through the re- 
sulting coherence minima. A profound change in the radial 
emission pattern is visible when one adopts the "reduced 
9" approach, which predicts significant radio emission up 
to higher distances depending on the relative orientation 
of observer and magnetic field. This is an important pre- 
requisite for the detection of ultra-high energy EAS with 
an array of affordable collecting area in combination with 
particle detector arrays such as KASCADE Grande or 
the Pierre Auger Observatory and will be verifiable by 
LOPES. 

The emitted total power in the coherent regime at 
low frequencies scales as the number of particles squared, 
which could therefore be probed directly by radio mea- 
surements of EAS, yielding information about the primary 
particle energy. 

We have not explicitly presented how variations of 
other parameters influence the radio emission, but most of 
the associated effects are fairly straight-forward to foresee: 
The emitted power scales linearly with B-field strength. 
The declination of the i?-ficld in Central Europe effectively 
decreases the value of B and introduces an asymmet- 
ric pattern to the radial dependence. An increase of the 
primary particle energy will boost the radiation because 
higher-energy showers will have their maximum closer to 
the observer. At the same time, the number of parti- 
cles increases linearly with primary particle energy and 
the power emitted at low frequencies increases as number 
of particles squared, which more than compensates the 
shrinking of the illuminated area on the ground. Inclined 
air showers will cause an asymmetric emission pattern and 
an attenuation of the emitted power because they reach 
their development at higher altitudes. A stronger curva- 
ture of the shower front will shift the interference minima 
to smaller radial distances and thus slightly decrease the 
effective size of the illuminated area on the ground. 

9.2. Comparison with experimental data 

A number of experiments have clearly established the pres- 
ence of radio emission from cosmic ray air showers in the 
past. A dependence of the polarisation of the emitted ra- 
diation on the earth's magnetic field direction was also 
confirmed by a number of experiments (e.g., Allan et al. 
1969), supporting the case for the geomagnetic emission 
mechanism. The actual strength of the emission, however, 
is still largely unknown at present state. The analysis of 



Allan (1971) led to a widely used formula summarising 
the presumed dependencies: 

= 20/.Vm-iMHz-i (^j^) 

X sin a cos ri exp ( ^ — - ) , (48) 

where the scale factor ro corresponds to (110±10) m at 
z/ = 55 MHz and for 77 < 35°. Later works (e.g., Sim 
1975; Prah 1971 and references therein), however, yielded 
values as low as 1-5 jiY m~^ MHz~^. A recent experiment 
in conjunction with the CASA/MIA array conducted by 
Green et al. (2003) was only able to place upper limits of 
= 31 — 34 fiV m~^ MHz~^ on the emission strength. 
Part of these discrepancies could be explained by un- 
certainties in the primary particle energy calibration at 
the time the experiments were made. A number of authors 
involved in the past works suspect the calibration of the 
radio measurements to be the major source of uncertainty 
(Atrashkevich et al. 1978). Additionally, the documenta- 
tion of the available data is not always totally precise re- 
garding the included energy ranges of primary particles, 
the selection of allowed zenith angles, the radial distance 
to the shower axis or the back-projection of the electric 
field vector in the plane normal to the shower axis and 
earth's magnetic field, which further complicates the is- 
sue. 

Extremely low values of of only 1 /xV m~^ MHz~^ or 
even lower are, however, disfavoured by the fact alone that 
air showers actually have been measured by experiments 
with only a few antennas (e.g., two in case of Prah 1971) 
with receivers of only a few MHz bandwidth in the early 
experiments. 

In this difficult situation, we choose to revert to the 
well documented data of Allan et al. (1970) as the ba- 
sis of our analysis. A comparison of these data with 
our predicted radial dependence of the emission is shown 
in Figure 17. While we clearly overpredict the emission 
strength in the centre, the general radial dependence fits 
relatively well. Regarding the spectral dependence, we 
make use of the Spencer (1969) data as presented, con- 
verted and complemented in Allan (1971) as well as the 
Prah (1971) data. These data sets, again, yield consider- 
ably lower values of Ci/, and we manually scale them up to 
make them consistent with the Allan et al. (1970) radial 
data. While the absolute values presented in Fig. 16 there- 
fore are somewhat arbitrary, the trend in the dependence 
actually does correspond to the spectral dependence that 
we predict near the shower core. 

All in all, we overpredict even the most optimistic past 
data by a factor 2, which is, however, not too surprising 
considering the very simplified integration over the shower 
evolution as a whole and the problems involved especially 
in the centre region. Additionally, the cutoff of the spa- 
tial integration as stated in Sec. 6.3 redistributes further 
emission to the centre region. 

We feel that having achieved a result which is consis- 
tent with past experimental data within a factor of "a few" 
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using such approximate descriptions of the shower char- 
acteristics and a mainly analytical approach incorporat- 
ing major approximations is a very encouraging outcome. 
In addition, our result further supports the geomagnetic 
emission mechanism as the dominant source of radio emis- 
sion from EAS. A huge improvement of our model will only 
be possible if we revert to elaborate computer simulations 
that use fewer approximations, more realistic particle dis- 
tributions (e.g., by interfacing our model to CORSIKA), 
and include additional effects such as the charge excess 
mechanism. Consequently, this will be the next step in 
our modeling efforts. 

Apart from further development towards a more so- 
phisticated model of radio emission from EAS, the most 
important aim for the near future therefore clearly is the 
obtainment of new, reliable data — as will be provided by 
LOPES, which should be able to easily measure the radio 
emission from a lO'^^ eV air shower as illustrated by the 
signal-to-noise levels overplotted in Figure 17. 

10. Conclusions 

We have analysed properties of radio emission from EAS 

in the scenario of coherent geosynchrotron emission. Our 
step-by-step analysis has helped to disentangle the co- 
herence effects arising from the different physical scales 
present in the air shower and to get a good feeling for 
the relative importance of these effects. While the spec- 
tral cutoff is directly governed by the longitudinal extent 
of the air shower, the radial dependence arises from the 
intrinsic beaming of the synchrotron radiation and its su- 
perposition over the shower evolution as a whole. 

The emitted radio power is of the expected order of 
magnitude, which is the strongest constraint we can make 
at the moment due to the large uncertainties associated 
with the available experimental data. Hence, in light of the 
data available to date, coherent geosynchrotron emission 
is able to explain the bulk of the radio emission from EAS. 

Our calculations show that LOPES should be able to 
easily detect the radio emission from a typical 10^^ eV 
air shower and will be a very useful tool for the study 
of EAS properties, especially the longitudinal structure of 
the particle distribution in the shower. 

In the future, we plan to use this model as a basis for 
the development of a sophisticated numerical computation 
including, among other aspects, 

— near-field effects (important for ultra-high energy EAS 
that develop their maximum near ground level) 

— vectorial integration for analysis of the polarisation of 
the radiation 

— Askaryan-type Cerenkov emission 

— and an interface to realistic air shower simulations such 
as CORSIKA. 
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The direction of the dominating emission component then 
changes as follows with (i?, <p): 



B X v{{), ip) 
I sina(^?, <^)| 



(55) 



These are all of the geometrical relations that are needed 
to execute the integration. 



Appendix: Geometry 

We adopt the instantaneous velocity vectors of the gen- 
erated particle pairs as radially pointing away from the 
centre of the sphere. For a particle at position (i?, ip) on 
the shell, its direction is therefore given by 

(cos r] sin 'd cos p + sinr] cos '& \ 
sin sin 1 , (49) 

sin r] sin i? cos ip — cos r] cos i? / 

whereas the direction of the B-field is given by 

sin rjB cos <Pb 

B= I sinryssint^B | ■ (50) 
— cos rjB 

Furthermore, the line-of-sight vector R from the particle 
to the observer is given by 

/ sin(r/ + j5io) \ 
R{^,^) = {Ro + K)\ \-Kv{^,^) (51) 

V - cos(?7 + do) J 

The direction of R is then calculated as 

and the pitch angle and angle to the line-of-sight corre- 
spond to 

cosa{i),Lp) = v{i3,(p) ■ B (53) 

cos6'(i?,v5) = v{{),ip)-h{{),ip). (54) 



